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A  i«<racf .Sliding  mode  control  is  examined  from  the  perspective  of  obtaining  stable  and  robust  tracking 
of  an  arbitrary  time-varying  reference  by  a  multi-input-output,  linear,  time-invariant  system  driven  by  a 
certain  class  of  bounded  errors,  nonlinearities  and  disturbances.  Most  existing  schemes  for  such  systems 
are  subsumed  by  the  one  presented  here.  The  results  are  developed  via  the  use  of  inverse  models,  and 
make  clear  the  constraints  imposed  by  the  finite  and  infinite  sero  structure  of  the  system.  In  particular, 
stable  and  robust  tracking  is  shown  to  be  obtained  by  the  scheme,  in  this  paper  if  and  only  if  the  system 
is  minimum  phase.  - -  ■  >  „•  *  „<.««  &(  €rC' ' 
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1.  Introduction 


Variable  structure  systems  with  sliding  modes  have  been  the  focus  of  a  growing  literature, 
see  for  example  [l]-[6]  and  references  therein.  The  results  in  this  paper  are  developed  in  the 
context  of  stable  and  robust  tracking  of  an  arbitrary  time- varying  reference  by  a  multivariable, 
linear,  time-invariant  system  driven  by  a  certain  class  of  errors,  nonlinearities  and  disturbances. 
Most  existing  sliding  mode  schemes  for  such  systems  are  subsumed  by  the  one  described  here. 
Our -treatment  is  built  around  the  use  of  inverse  models,  and  offers  some  new  insights  into 
the  constraints  imposed  by  the  finite  and  infinite  zero  structure  of  the  system.  It  is  shown 
that  stable  and  robust  tracking  is  obtained  by  this  sliding  scheme  if  and  only  if  the  system  is 
minimum  phase.  The  blending  of  state-space  and  polynomial  matrix  descriptions  allows  our 
results  to  be  derived  cleanly  and  generally. 


Stable,  Robust  Tracking.  Our  starting  point  is  the  system 


x(t)  =  Az(t)  +  Bu(t)  +  v(x,  t ) 


*  Work  of  the  first  author  was  supported  by  the  Army  Research  Office  under  Grant  DAAG-29-84-K- 
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where  z  is  the  n-dimensional  state  vector  and  u,  y  are  m-dimensiona)  control  input  and  tracking 
output  vectors  respectively.  The  vector  v  represents  model  errors,  nonlinearities,  and  distur¬ 
bances.  (The  time  argument  t  will  be  dropped  for  notational  simplicity  wherever  it  is  not 
needed  for  clarity  or  emphasis.)  One  can  obviously  not  ask  for  more  independently  tracking 
outputs  than  there  are  control  inputs.  On  the  other  hand,  if  the  number  of  desired  tracking 
outputs  is  fewer  than  the  number  of  control  inputs,  then  there  are  additional  degrees  of  freedom 
that  may  be  exploited,  and  some  remarks  on  the  possibilities  are  made  in  the  Conclusion  of 
this  paper. 

The  control  objective  is  to  pick  u  as  a  function  of  the  measured  z  so  as  to  have 

«(0  =  v(0  -  y^(0  (2) 

go  to  rero,  where  yd  is  the  desired  output,  i.e.,  the  tracked  reference,  and  is  assumed  to  be 
sufficiently  differentiable  —  the  precise  degree  of  differentiability  will  become  explicit  later. 
We  also  require  that  the  dynamics  of  z  be  independent  of  v  —  which  gives  rise  to  the  label 
“robust”  —  and  that  when  yd  —  0  identically  then  z  goes  asymptotically  to  rero  —  hence  the 
label  “stable”. 

Assumptions.  Assume  that  the  system  (l)  is  controllable  by  u  and  observable  from  y. 
Also  impose  the  standard  matching  restriction,  (2], [4],  on  v,  namely 

v{x,t)  -  Bw(x,t)  (3) 

for  some  w  that  has  a  known  bound  in  some  norm.  It  is  these  assumptions  on  v  that  end  up 
enabling  robust  tracking.  Finally,  assume  that  the  transfer  function  matrix 

G(A)  =  C(\I  -  A)~l B  (4) 

of  the  system  (1)  is  invertible  (as  a  rational  matrix).  This  assumption  is  clearly  necessary  if 
tracking  of  arbitrary  references  is  desired. 

The  invertibility  assumption  on  the  transfer  function  implies  that  it  has  as  many  zeros  as 
poles,  if  seros  at  infinity  are  included  and  if  multiplicities  are  properly  accounted  for,  [7].  The 
controllability  and  observability  of  the  model  imply  that  this  number  (of  poles  and  of  seros) 
is  n,  the  order  of  the  system  (1).  The  finite  seros  are  (see  [8,  Sec.  6.5])  precisely  the  roots  of 
det(A/-  A).detC(A).  If  t  denotes  the  relative  order  of  detG(A),  i.e.  its  denominator  degree 
minus  its  numerator  degree,  it  follows  that  there  must  be  n  -  t  finite  seros.  Hence  tht  total 
number  of  infinite  seros  is  t. 
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Sketch  of  Main  Results.  First  consider  the  special  case  where  CB  is  invertible.  Suppose 
the  control  u  is  chosen  to  bring  the  system  in  finite  time  from  any  initial  state  to  the  (n  -  tri¬ 
dimensional  surface  given  by  s(x, t)  =  Cx -y<j  =  y-y<<  =  e  =  0  and  to  then  keep  the  trajectory 
confined  to  —  or  "sliding”  on  —  this  surface  by  switching  appropriately,  so  that  e  =  0  from 
then  on.  That  this  can  always  be  done  under  the  above  assumptions  is  well  established  in  the 
literature  on  sliding  modes,  see  [l]-{6],  so  tracking  is  indeed  achieved  in  this  case. 

To  determine  the  input  and  state  trajectories  when  the  system  is  sliding,  one  can  use  the 
method  of  equivalent  control  [l],(3].  This  effectively  sets  ds/dt  =  0  for  motion  on  the  sliding 
surface  and  solves  for  uM,  which  is  the  Filippov  average  of  the  switched  control,  [6]: 

ds/dt  =  y  -  i/d  =  Cx  -  Vd  =  CAx  +  CB(uM  +  u>)  -  yd  =  0  (5) 

from  which,  since  CB  is  assumed  invertible, 

=  -w  -  (CB)~l(CAx  -  i/d)  (6) 

The  system  behavior  in  the  sliding  regime  is  therefore  given  by 

*  =  Ax  +  £(u«t  +  u>) 

*  (/  -  B{CB)~lC)Ax  +  B{CB )-xfe,  (7) 

and  is  unaffected  by  w.  Note  that  we  have  not  guaranteed  robustness  to  w  during  the  interval 
when  the  system  trajectory  is  moving  towards  the  sliding  surface  but  not  yet  sliding.  This 
problem  can  be  overcome  by  ensuring  that  the  system  starts  in  the  sliding  mode,  and  this  can 
be  achieved  by  choosing  1/4(0)  =  Cz(0). 

Though  we  have  managed  to  achieve  robust  tracking  in  the  sliding  mode,  the  dynamics  in 
(7)  is  evidently  fixed,  and  may  be  unsatisfactory.  The  characteristic  polynomial  of  the  system 
in  (7)  is  easily  seen  to  be 

det(A I-A  +  B{CB)-lCA]  =  det(C£)"1.Am.det(A/-  A).detG(A)  (8) 

Note  from  this  that  the  matrix  in  (7)  has  m  eigenvalues  at  0,  corresponding  to  the  fact  that  the 
trajectory  is  confined  to  an  (n  -  m)-dimensional  surface.  The  remaining  eigenvalues  are  seen 
from  (8)  to  occur  at  the  n  -  m  finite  seros  of  the  transfer  function  G(A).  This  may  have  been 
anticipated  by  considering  the  case  where  1/4  =  0  for  all  time:  evidently  the  equivalent  control 
must  in  this  case  be  exciting  the  system  at  the  frequencies  of  its  seros  in  order  to  keep  y  =  0 
when  “««  +  o-  (This  observation  has  also  been  made  in  [9].)  Hence  stable  tracking  is  obtained 
with  this  scheme  if  and  only  if  the  system  is  minimum  phase. 


When  CB  ia  not  invertible,  uM  cannot  be  determined  by  simply  taking  the  first  derivative 
of  y  as  in  (5).  A  more  elaborate  procedure,  familiar  from  the  area  of  system  inversion,  (10)- 
[12],  is  now  required.  The  procedure  in  effect  takes  appropriate  linear  combinations  of  further 
derivatives  of  the  components  of  y  to  compute  u(f .  However,  the  equivalent  control  method 
requires  that  all  of  this  procedure  be  folded  into  the  simple  act  of  computing  d$/dt.  This 
implies  that  the  sliding  surface  above  has  to  be  modified,  and  we  shall  show  that  it  goes  from 
8  =  e  =  Cx  -  yd  =  0  to  s  =  P(d/dt)e  =  Lx-  P(d/di)yd  =  0  for  some  appropriately  chosen 
polynomial  matrix  P( A)  and  corresponding  coefficient  matrix  L. 

It  will  turn  out  under  our  assumptions  that  LB  is  invertible,  so  the  earlier  analysis  can 
be  followed,  with  C  replaced  by  L  and  with  other  obvious  changes.  Now  the  tracking  error  e 
in  the  sliding  regime  has  dynamics  that  is  governed  by  the  equation 

P(d/dt)e  =  0  (9) 

The  associated  characteristic  polynomial  det P(\)  will  be  shown  to  have  degree  t-m  and  to 
be  arbitrarily  assign  able.  We  shall  also  demonstrate  that  the  dynamics  of  x  is  unaffected  by 
xv  in  the  sliding  mode.  The  eigenvalues  governing  the  state  dynamics  will  be  seen  to  comprise 
precisely  the  t—m  roots  of  det  P( A),  the  n  —  t  finite  zeros  of  G( A),  and  m  eigenvalues  at  0  that 
correspond  as  before  to  the  confinement  of  the  trajectories  to  the  (n  -  m)-dimensional  sliding 
surface.  Again,  the  appearance  of  the  finite  zeros  of  the  system  may  have  been  anticipated  from 
the  fact  that,  with  yd  =  0  identically,  the  output  y  can  be  made  to  go  to  0  arbitrarily  fast  in 
this  scheme. 

Outline  of  Paper.  The  remainder  of  the  paper  fills  out  and  proves  the  above  results. 
Section  2  contains  the  main  results.  Inverse  models  are  introduced  in  a  way  that  is  convenient 
for  our  purposes,  and  their  role  in  defining  sliding  surfaces  and  the  associated  sliding  dynamics 
is  established.  This  section  also  briefly  considers  how  to  pick  control  laws  that  cause  trajectories 
to  reach  a  sliding  surface  in  finite  time.  Section  3  presents  some  illustrative  examples.  The 
Conclusion  of  the  paper  mentions  possible  extensions.  Certain  proofs  whose  details  are  not 
central  to  the  flow  of  the  paper  are  gathered  in  Appendices  A1-A3. 

2.  Sliding  Modes  Defined  Through  Inverse  Models 

Inverse  Models.  The  literature  on  system  inversion,  see  [10j-[12]  for  instance,  shows  how 
one  may  compute  (via  the  structure  algorithm  in  [10]  for  example)  an  invertible  polynomial 
matrix  Q(A)  such  that 

awcmI,...-  7  (>°) 

In  the  language  of  [13],  Q(A)  is  an  ‘identity  interactor”  for  the  system  (l)  or  the  transfer 
function  (4).  The  reason  that  Q( A)  appears  explicitly  or  implicitly  in  all  approaches  to  system 


inversion  is  that  the  operator  Q(d/dt),  acting  on  the  system  output,  takes  linear  combinations 
of  derivatives  of  the  components  of  this  output  vector  to  produce  a  signal  that  directly  contains 
the  input.  This  is  suggested  by  (10)  but  will  be  verified  more  precisely  by  (14)  below. 

It  can  be  seen  from  (10)  that  the  relative  order  of  detQ(A).detG(A)  is  0.  Since  the 
relative  order  of  G(A)  is  l,  it  follows  that  the  degree  of  the  polynomial  detQ(A)  is  also  t.  It  is 
further  true  that  the  detailed  polar  structure  of  Q(A)  at  A  =  oo,  in  the  sense  of  [7],  is  identical 
to  the  zero  structure  of  G(A)  at  A  =  oo  or  the  polar  structure  of  G-1(A)  at  A  =  oo.  In  particular, 
the  highest  power  of  A  in  Q( A)  equals  its  highest  power  in  G-l(A). 

The  Q(A)  in  (10)  is  not  unique,  however.  It  is  not  hard  to  show  that  if  0(A)  is  one 
particular  solution  of  (10),  then  a  polynomial  matrix  0(A)  is  a  solution  if  and  only  if 

0(A)  =  0(A) +  JNT(A)  (11) 

where  N(X)  is  a  polynomial  matrix  such  that  N( X)Q~1(X)  is  strictly  proper  (i.e.  vanishes  at 
A  =  oo).  A  necessary  condition  on  N( A)  for  this  strict  properness  is  that  its  column  degrees  (i.e. 
the  highest  degree  among  the  entries  in  each  of  its  columns)  be  less  than  the  corresponding 
column  degrees  of  0(A).  It  follows  from  (11)  that  all  solutions  0(A)  have  the  same  column 
degrees.  This  simple  condition  for  strict  properness  is  also  sufficient  if  and  only  if  0(A)  is 
"column  reduced”,  i.e.  if  and  only  if  the  degree  l  of  det0(A)  equals  the  sum  of  the  column 
degrees  of  0(A),  see  [8,  Sec.  6.3]. 

We  now  examine  more  precisely  the  role  of  0(A)  in  defining  an  inverse  model  for  the 
given  system  (1),  (3).  It  is  shown  in  Appendix  A1  that 

J-0(A)G(A)  =  ff(A/-A)-1B  (12) 

for  some  constant  matrix  H.  The  controllability  of  the  system  guarantees  that  (XI  —  A)~lB 
has  no  constant  vectors  in  its  left  nullspace,  so  H  can  be  found  as  the  unique  constant  solution 
of  (12)  for  a  given  0(A)*  (Another  route  to  computation  of  H  is  provided  by  the  polynomial 
matrix  division  described  in  Appendix  Al,  while  a  third  approach  is  contained  in  (23b)  below.) 
On  application  of  the  well-known  matrix  inversion  lemma  to  (12),  one  obtains  the  following 
expression  for  the  transfer  function  of  the  inverse  system 

G-l(A)  =  [£T(A/-  A  -  BH)~lB  +  J]Q(A)  (13) 

Note  that  the  seros  of  G(A)  at  finite  values  of  A  give  rise  to  poles  of  G-l(A)  at  the  same 
locations.  Since  Q( A)  has  no  finite  poles,  (13)  shows  that  these  poles  must  be  eigenvalues  of 
A  +  BH.  In  other  words,  the  eigenvalue*  of  A  +  BH  mutt  include  all  the  finite  zero*  of  G(A). 


Given  the  expression  in  (13)  for  the  inverse  of  the  system  transfer  function,  it  is  perhaps 
not  surprising  that  the  input  u  +  w  of  the  system  (l),  (3)  can  in  fact  be  determined  from  the 
output  y  via  the  inverse  model 

i  =  {A  +  BH)x  +  BQ{d/dt)v  (14o) 

u  + to  =  Hx  +  Q{d/dt)y  (146) 

where  the  state  vector  x  is  the  same  as  that  in  (1).  A  proof  of  this  is  given  in  Appendix  A3.  The 
inverse  model  (14)  takes  as  its  input  the  signal  Q(d/dt)y,  which  constitutes  linear  combinations 
of  derivatives  of  components  of  the  output  of  the  original  system  (1).  The  output  of  the  inverse 
model  (14)  is  the  input  u  +  w  of  the  original  system  (1),  (3).  This  confirms  the  claim  following 
(10)  regarding  the  significance  of  Q(A). 

The  Sliding  Surface.  The  inverse  model  (14)  shows  that  u  +  w  can  be  computed  from 
Q(d/dt)y.  Recall  from  the  sketch  in  the  Introduction  that  the  equivalent  control  method  com¬ 
putes  by  setting  dafdt  —  0.  This  suggests  using  a  sliding  surface  with  the  property  that 
da/dt  =  Q(d/dt)e,  which  can  be  achieved  by  picking 

a  =  P(d/dt)e  =  0  (15) 

where  P( A)  is  a  polynomial  matrix  that  satisfies 

AP(A)  =  Q(A)  (16) 

It  is  obvious  from  (16)  that  such  a  P( A)  exists  if  and  only  if 

9(0)  =  0  (17) 

Since  the  degree  of  detQ(A)  is  l,  that  of  detP(A)  is  seen  from  (16)  to  be  l-m.  The  number 
of  degrees  of  freedom  available  in  picking  P( A)  and  hence  the  sliding  surface  will  be  shown  in 
the  next  subsection  to  be  m[t  -  m). 

With  s  chosen  as  in  (15),  (16), 

da/dt  =  ( d/dt)P(d/dt)e  =  Q(d/dt)e  =  Q(d/dt)y  -  Q{d/dt) yd  (18) 

as  desired.  Now  setting  da/dt  =  0  in  the  sliding  mode,  as  prescribed  by  the  equivalent  control 
method,  and  using  the  inverse  model  (14b)  to  substitute  for  Q{d/dt)y,  we  find  that 

s=-w  +  Hx  +  Q(d/dt)yd  (19) 

With  this  control  input  in  (1),  (3),  the  state  dynamics  in  the  sliding  mode  is  governed  by 

x  =  (A  +  BH)x  +  BQ(d/dt)yd  (20) 
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and  the  output  tracking  error  satisfies  (15),  (9).  Once  again,  therefore,  we  obtain  robust 
tracking  after  sliding  is  established. 

Before  proceeding  further,  we  need  to  show  that  (15),  (16)  does  actually  define  a  surface 
in  state  space.  Appendix  AS  demonstrates  that 

P[d/dt)y  =  P(d/dt)Cx  =  Lx  (21) 

for  a  constant  matrix  L  that  is  uniquely  determined  by  P( A),  so 

s(x,  t)  =  P[d/dt)y  -  P(d/dt)yd  =  Lx-  P(d/dt)yd  (22) 

Hence  (15)  does  indeed  constitute  a  surface  in  state  space.  To  compute  L  from  (21),  one 
evaluates  P(dJdt)Cx,  with  repeated  elimination  of  x  using  (la)  —  this  is  what  the  polynomial 
matrix  division  in  Appendix  Al,  used  in  Appendix  A3,  effectively  does.  The  derivation  in 
Appendix  A3  guarantees  that  u  +  vj  and  its  derivatives  will  not  appear  in  the  result,  so  the 
final  result  is  of  the  form  Lx. 

The  matrices  H  and  L  are  seen  from  Appendices  Al  and  A3  to  be  related  by 


LB  —  l 

(23a) 

i 

II 

hj 

(236) 

so  £  is  a  left  inverse  of  B.  It  should  not  be  inferred  that  every  choice  of  left  inverse  L  leads 
to  a  sliding  surface.  Indeed,  we  have  already  claimed  above  (and  will  show  below)  that  there 
are  only  m(t  —  m )  degrees  of  freedom  in  picking  the  sliding  surface,  while  there  are  m(n  —  m) 
degrees  of  freedom  in  picking  left  inverses  of  B.  Once  L  is  found  from  (21)  or  an  equivalent 
computation,  then  H  may  be  found  from  (23b),  as  mentbned  earlier. 

With  «  identified  as  in  (22),  the  development  in  (5)-(7)  can  be  followed,  with  C  and  yd 
replaced  by  L  and  P[d/dt)yd  respectively.  Using  the  relations  in  (23),  the  result  is  again  (20), 
as  expected,  which  is  the  generalization  of  (7).  Note  that  in  order  to  ensure  that  the  system 
starts  directly  in  the  sliding  mode,  one  must  pick  yd  and  its  derivatives  such  that  s(x(0),0)  =  0. 

Dynamics  in  the  Sliding  Mode.  Without  loss  of  generality,  the  tracked  reference  yd 
can  and  will  be  taken  to  be  0  when  dealing  in  this  subsection  with  stability  and  dynamics.  The 
state  dynamics  (20)  in  the  sliding  regime  is  then  undriven.  Note  from  (15)  that  the  dynamics 
of  the  tracking  error  in  the  sliding  mode  is  governed  by  the  matrix  P(A)  and  the  associated 
characteristic  polynomial  det  P( A),  whose  degree  we  have  already  identified  as  t  -  m.  Taking 
into  account  the  parametrization  of  <J(A)  in  (11)  and  the  relation  between  P( A)  and  Q(X)  in 


(16),  it  is  teen  that  a  complete  parametrizaticn  of  P(A)  in  terms  of  a  particular  solution  P(A) 
is  given  by 

P(A)  =  P(A)  +  M(A)  (24) 

where  M( A)  is  a  polynomial  matrix  such  that  M(\)P~l(\)  is  strictly  proper.  Since  the  column 
degrees  of  M( A)  are  less  than  the  corresponding  degrees  of  P(A),  it  follows  that  the  column 
degrees  of  P( A)  are  invariant.  Hence  the  orders  of  the  highest  derivatives  of  the  components  of 
e  required  to  form  the  sliding  surface  are  invariant. 

The  question  now  is  what  freedom  there  is  in  assigning  the  tracking  error  dynamics. 


Combining  results  in  [8,  Sec.  6.3,  Sec.  7.2]  with  the  parametrization  in  (24)  shows  that  (24) 
affords  us  m(t- m)  degrees  of  freedom,  and  that  we  can  in  particular  use  this  freedom  to  make 
(in  many  ways,  if  m  >  1)  the  polynomial  detP(A)  equal  to  any  polynomial  of  degree  t-  m, 
up  to  a  nonzero  constant  multiple.  Arbitrarily  fast  decay  of  the  output  tracking  error  may 
therefore  be  obtained  (assuming,  as  in  our  problem  formulation,  that  the  control  input  is  not 
constrained). 

We  turn  from  examination  of  the  output  tracking  error  dynamics  to  that  of  the  system 
state  in  the  sliding  regime.  Since  y  cannot  display  dynamic  behavior  that  is  not  present  in  z, 
it  must  be  true  that  the  eigenvalues  of  A  +  BH  in  (20)  include  all  the  l  -  m  roots  of  det  P( A). 
Also,  we  have  already  argued  that  the  matrix  A  +  BH  must  contain  all  the  n  —  l  finite  zeros  of 
the  given  system  (1)  among  its  eigenvalues.  Finally,  A  +  BH  must  have  m  eigenvalues  at  the 
origin,  corresponding  to  the  confinement  of  the  trajectories  to  the  (n-m)-dimensional  subspace 
Lx  =  0.  All  of  this  therefore  accounts  for  [£  —  m]  +  [n  - £]  +  m  —  n  eigenvalues  of  A  +  BH,  and 
there  are  no  more  to  account  for!  This  can  be  explicitly  and  directly  shown  by  verifying  that 


(25) 


det(A/-  A  -  BH)  =  detP(A).Am.det(A/-  A).detG(A) 

using  the  various  relationships  already  established.  In  conclusion: 

The  dynamic  behavior  of  the  system  on  the  gliding  turf  ace  it  governed  by  n— l  eigenvalues 


that  are  located  at  the  finite  zerot  of  the  tyttem  and  a  further  l—m  eigenvaluet  that  are 
freely  assignable,  subject  to  being  a  self-conjugate  set. 

It  follows  that  stable  tracking  is  possible  with  this  scheme  if  and  only  if  the  system  is  minimum 
phase. 

Obtaining  Sliding.  The  development  so  far  has  focused  on  what  happens  once  the  sliding 
mode  is  attained.  In  order  to  see  what  sort  of  nonlinear  control  law  is  needed  to  reach  and  stay 
on  the  sliding  surface,  note  from  (14b)  and  (18)  that 


ds/dt  =  Q{d/dt)y  -  Q{d/dt)yd  =  u  +  w  -  Hx  -  Q[d/dt)yd 
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(26) 


Now  pick  the  control  law  to  be 


u  =/(•)  + Hx  +  Q(d/dt)  yd  (27) 

so  that 

da/dt=  f(a)  +  w  (28) 

The  function  /(«)  can  be  picked  so  as  to  guarantee  convergence  of  a  in  (28)  to  0  in  finite  time 
in  the  presence  of  a  bounded  w.  Many  ways  for  doing  this,  using  switching  control  laws,  are 
discussed  in  the  sliding  mode  literature,  and  so  will  not  be  considered  here.  It  should  be  noted 
that  our  formulation  has  not  considered  the  bounds  on  control  that  will  exist  in  any  realistic 
situation,  and  these  bounds  must  be  explicitly  reckoned  with  in  any  actual  design.  The  effect 
of  the  choice  of  Q( A)  (and  hence  H),  of  Vd,  and  of  f(a)  on  the  control  u  must  be  assessed  from 
(27). 

One  final  point  is  worth  making.  The  sliding  surface  specified  by  a  =  0  can  also  be 
specified  by  the  condition  Ta  =  0  where  T  is  any  nonsingular  m  x  m  matrix.  Working  with  the 
sliding  surface  specified  in  this  form  does  not  change  any  of  the  expressions  for  motion  while 
sliding,  though  it  will  change  the  expressions  that  govern  the  approach  to  the  sliding  surface. 

3.  Examples 

Example  1.  If  CB  is  invertible,  the  system  (1)  has  1=  m  infinite  zeros  and  n  -  m  finite 
ones.  The  set  of  Q(X)  satisfying  (10),  (11)  is  given  precisely  by  Q(A)  =  (CB)~1(XI+  N),  where 
N  is  an  arbitrary  constant  matrix.  Imposing  the  additional  constraint  in  (17)  that  Q(0)  =  0 
yields  Q( A)  =  (CB)'1  A  as  the  only  solution.  The  corresponding  P( A)  is  simply  (CB)'1.  It 
follows  that  the  sliding  surface  can  be  defined  via  (CB)~l(Cx  -  yd)  =  0  or  equivalently  via 
Cx  -  yd  =  0,  which  is  the  same  surface  as  the  one  considered  in  the  Introduction.  The  results 
of  the  Introduction  are  recovered  on  applying  (l5)-(20).  •• 

Example  2.  Suppose  the  model  (l)  has  order  n  =  2m,  and  that  its  coefficient  matrices 
are  given  by 

A=  B  =  ( j-i )  >  C  =  (Im  0)  (29) 

This  model  represents  the  prototype  spring-inertia  system,  with  J,  K  and  D  being  inertia, 
spring  and  damping  coefficient  matrices  respectively.  Such  models  are  found  in  descriptions 
of  robot  manipulator  dynamics  [5], [6],  in  the  swing  equations  that  describe  electromechanical 
oscillations  of  power  systems,  and  so  on.  The  first  block  of  state  variables  comprises  position 
variables  6 ,  which  constitute  the  outputs  in  this  example,  the  second  block  comprises  velocity 
variables  v,  and  the  inputs  are  generalized  forces. 
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For  this  model, 


<7(A)  =  (AV  +  £>A  +  if)-1 


(27) 


All  n  =  2m  zeros  of  this  system  are  at  oo,  so  t  =  2m.  The  set  of  Q(A)  in  (10),  (11)  is  precisely 
given  by  Q[ A)  =  J(\iI+M\  +  N),  where  M  and  N  are  arbitrary  constant  matrices.  Imposing 
the  constraint  (17)  yields  N  =  0.  We  then  have  P( A)  =  J(A/+  M).  The  roots  of  detP(A)  are 
evidently  the  m  eigenvalues  of  M,  and  can  be  assigned  to  any  self-conjugate  set  of  values.  The 
corresponding  sliding  surface  is  given  by  (w  -  6d)  +  M{8  -  Sd)  =  0.  •• 

Example  3.  Let 


/  0  1  0  (A  (0  0\ 

-1  -2  0  0  1  0 

0  0  — 1  0  ’  15  —  0  1 

^ooooj  [o  lj 


(28  a) 


C=(22  1  1  5) 

C(A'= ({a+ i/(*/+i))  (29> 

Using  the  methods  described  in  [7],  this  transfer  function  matrix  can  be  seen  to  have:  1  first- 
order  pole  at  A  =  0;  1  first-order  and  1  second-order  pole  at  A  =  -1;  1  first-order  zero  at  A  =  -2; 
and  1  first-order  and  1  second-order  zero  at  A  =  00.  The  total  number  of  infinite  zeros  t  is  thus 
3.  One  can  show  that  the  general  form  for  Q(A)  is 

nfW—  (~ A3  +  »iA  +  n3  A2  +  (1  -  nj)A  +  n6  ^  ,,n> 

^"lA’  +  njA+n,  -A2-nsA  +  n«  )  ^ 

for  arbitrary  n,-.  Imposing  the  condition  in  (17)  and  computing  P(A)  according  to  (16)  yields 

pW  =  ('a++«"‘  <31> 

The  determinant  of  this  matrix  is  -(A  +  ns),  which  is  of  degree  1  and  arbitrarily  assignable,  as 
expected.  The  corresponding  sliding  surface,  taking  nt  =  0  for  instance,  can  then  be  computed 
from  (21),  (22)  to  be 


Conclusion 

Most  existing  approaches  to  sliding  mode  control  have  not  considered  the  tracking  prob¬ 
lem  posed  here.  The  present  perspective  is  more  useful  in  many  situations,  and  includes  the 
main  results  developed  via  other  approaches.  If  no  tracking  outputs  are  specified  a  priori,  then 
one  is  free  to  choose  C  to  obtain  desirable  system  dynamics.  This  is  the  form  in  which  the 
sliding  mode  control  problem  is  more  commonly  considered  in  the  literature.  The  results  in 
this  paper  show  that  C  should  be  picked  to  place  any  finite  transmission  zeros  at  desirable 
locations.  Rosenbrock’s  theorem  [14,  Ch.5,  Thm.4.1]  indicates  the  freedom  available  in  doing 
this.  If  some,  but  fewer  than  m,  tracking  outputs  are  designated,  one  has  the  nontrivial  problem 
of  augmenting  the  output  set  to  obtain  a  C  that  leads  to  zeros  at  desirable  locations.  Some 
discussion  of  such  problems  is  contained  in  [15].  Note  that  any  zeros  of  the  original  system 
before  augmentation  will  remain  as  zeros  of  the  augmented  system. 

One  can  examine  extensions  of  the  results  here  to  the  nonlinear  case,  building  on  the 
notions  of  left  invertibility  for  such  systems,  [16].  Exploration  in  this  direction  may  be  found 
in  [17]. 

Appendix  Al.  Note  from  (10)  that 

Q(A)C(AI-A)~1B  =  I+B(A)  (Al.l) 

for  some  strictly  proper  R(X).  Carrying  out  division  [8,  Sec.  6.3]  of  Q{X)C  by  XI  -  A  yields 

Q(X)C  =  L(A)(A/  -  A)  -  H  (A1.2) 

where  L(A)  is  a  polynomial  matrix  and  H  is  a  constant  matrix,  both  uniquely  determined.  Now 
(Al.l)  becomes 

[L(A)-#(A/-  A)-I]B  =  /  +  R(A)  (A1.3) 

liquating  strictly  proper  parts  on  both  sides  of  (A1.3)  and  substituting  for  f2(A)  in  (Al.l) 
establishes  (12).  Equating  polynomial  parts  shows  that  L( A)  is  a  polynomial  left  inverse  of  B: 

L(X)B  =  I  (A1.4) 

Denoting  L( 0)  by  L,  (23a)  follows  from  (A1.4). 

Appendix  A2.  Using  (J),  (3),  (A1.2)  and  (A1.4),  we  see  that 

Q{d/dt)y  =  Q[d/dt)Cx  =  L{d/di)\I(d/dt)  -  A]z  -  Hx 

=  L(d/dt)B(u  +  tu)  -  Hx  =  (u  +  w)  —  Hx 


or  u  +  w  =  Hx  +  Q[d/dt)y ,  which  is  (14b),  and  which  yields  (14a)  oo  substitution  in  (1). 
Appendix  A3.  From  (16)  and  (A1.2),  and  recalling  that  L  =  1.(0),  it  follows  that 

P(A)C  =  A"1  [1(A)  -  L](A  J  -  A)  +  L  -  X~ 1  [LA  +  H )  (A3.1) 

Now  (17)  and  (A1.2)  show  that 

LA+H  =  0  (A3.2) 

(which  is  (23b))  so  the  last  term  in  (A3.1)  drops  out.  Also  note  that  the  term  A-1  [L(A)  —  L\  is 
polynomial.  From  (A1.4)  we  see  that 

(L(A)  -  L]B  =  0  (A3.3) 

Putting  together  (1)  and  (A3.l)-(A3.3)  yields  (21). 
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